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In this paper we operate under the assumption that no tensors from inflation will be measured in 
the future by the dedicated experiments and argue that, while for single-field slow-roll models of 
inflation the running of the spectral index will be hard to be detected, in multi-held models the 
running can be large due to its strong correlation with non-Gaussianity. A detection of the running 
might therefore be related to the presence of more than one active scalar degree of freedom during 
inhation. 


Motivations .— Inflation [1] provides a mechanism to ex¬ 
plain the initial conditions for the Large Scale Structure 
(LSS) and for the Cosmic Microwave Background (CMB) 
anisotropy. The seeds for the density scalar perturbations 
are created from the quantum fluctuations “redshifted” 
out of the Hubble radius during the early period of su¬ 
perluminal expansion of the universe. At the same time, 
the generation of tensor (gravity-wave) fluctuations is a 
generic prediction of an accelerated de Sitter expansion 
of the universe [2]. 

To characterize the scalar and tensor perturbations 
several observables are introduced: the spectral index (or 
tilt) ns of the scalar perturbations, currently measured 
to be Us = 0.9655 ± 0.0062 at 68% CL; the running with 
the wavenumber k of the spectral index as = dn^/dlnfe, 
currently consistent with zero, as = —0.00571 ± 0.0071 
at 68% CL (assuming no tensors); and the level of lo¬ 
cal non-Gaussianities parametrized by /nl = 0.8 ± 5.0 
at 68% CL and 5 nl = (-9.0 ± 7.7) x 10^ at 68% CL. 
All these figures are quoted from the last Planck releases 
[3, 4]. Finally, the tensor-to-scalar ratio is bounded to be 
r < 0.12 at 95% CL [5]. 

Many next-generation satellite missions, see for in¬ 
stance Refs. [6-8], are dedicated to measure the polar¬ 
ization of the CMB anisotropies. On large angular scales 
the B-mode polarization of the CMB carries the imprint 
of primordial gravitational waves, and its precise mea¬ 
surement would provide a powerful probe of the epoch 
of inflation. The goal of these missions is to achieve a 
measurement of r down to 0(10“^). 

While a detection of gravity waves from inflation will 
be a major milestone for the inflationary paradigm, on 
the other hand one needs to remember that going from 
r = 10“^ down to r = 10“^ means testing the energy 
scale of inflation Ainf only by a factor (100)^/^ ~ 3.1 
better than current bounds (r scales like 

One should therefore envisage the situation in which 
Nature has not be so kind to us to set the scale of the 
energy of inflation within a factor of three away from the 
current limits. In this paper we will therefore operate 
under the working assumption that no tensors will be 


measured in the future by the dedicated experiments. 

While we hope this hypothesis will prove to be wrong 
in the future, we believe it is reasonable to ask how else 
we could learn something more about the inflationary 
perturbations. 

In this paper we turn our attention to the scale de¬ 
pendence of the tilt of the scalar perturbations. Our hy¬ 
pothesis on the tensor modes leads to a hardly detectable 
running within single-field slow-roll models. On the other 
hand, a negative result on tensors might indicate a low 
level of the Hubble rate, leaving open the possibility that 
the scalar perturbations are due to a light scalar field 
other than the inflaton field. In such a case, our findings 
indicate that a large running may be achieved thanks to 
an interesting correlation between the running and the 
level of non-Gaussianity. Such strong correlation exists 
only in multi-field models of inflation, while it is miss¬ 
ing in single-field slow-roll models. If no tensors will 
be observed, a measurement of a sizable running of the 
spectral index might therefore imply a sizable amount 
of non-Gaussianity (and viceversa) and the presence of 
more than one active scalar degree of freedom during in¬ 
flation. 

Our considerations are particularly relevant for the 
goals of next missions designed to measure the three 
dimensional structure of the universe, such as the ESA 
mission EUCLID [9] and the NASA SPHEREx mission 
[10], a proposed all-sky spectroscopic survey satellite. 
Eor instance, the combination of the galaxy power spec¬ 
trum and bispectrum leads to the forecasted 68% CL 
errors a{ns — 1) = 2.2 x 10“^, (j{as) = 6.5 x 10“^ and 
^■(/nl) = 0.2 (for the local case) in the case of SPHEREx. 
This is an improvement with respect to EUCLID fore¬ 
casts of a factor ~ 2 in Ofs and a factor ~ 20 for the local 
/nl (or ~ 6 with the galaxy power spectrum only). For 
the other parameters SPHEREx will deliver constraints 
comparable to those expected from the EUCLID spec¬ 
troscopic survey. 

As EUCLID and (possibly) SPHEREx will strongly 
narrow down the allowed parameter space in the 
(os )/nl (-plane, it is therefore extremely interesting and 
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timely to understand what are the theoretical expec¬ 
tations in terms of and its connections to non- 
Gaussianity. 

Single-field slow-roll inflation .— Let us first consider the 
case in which the cosmological perturbations are gener¬ 
ated by the same scalar field </> driving inflation. The 
power spectrum of scalar perturbations is given in such 
a case by 
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2Mp2e* 
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where the sub-index * indicates that quantities have to be 
computed at Hubble crossing, iJ* is the Hubble rate dur¬ 
ing inflation, Mp ~ 2.4 x 10^® GeV is the reduced Planck 
mass, and e* = —is one of the slow-roll parame¬ 
ters (dots indicate derivatives with respect to time). 

Let us now assume that no tensors will be observed 
in the future and see what are the consequences we can 
draw. Measuring no tensors implies that r = 16e* ^ 
10“® or e* ^ 10“^. Irrespectively of the form of the 
inflaton potential, this means that ^ Mp. To see 
this, take the cosmological range of scales to span four 
decades, corresponding to In /c ~ 10. This corresponds 
to ~ 10 e-folds of inflation. In first approximation in 
slow-roll inflation e* has negligible variation over one e- 
fold and in typical models it has only small variation 
over the ~ 10 e-folds. Taking that to be the case, one 
finds (this is nothing else that the reverse of the so- 
called Lyth bound [11]) 0 <C 0.4 Mp, where we have used 
(primes indicate derivative with respect to the inflaton 
field) iV'(0) = ±(l/-\/^Mp) and N is the number of 
e-folds till the end of inflation. A more detailed compu¬ 
tation leads to sub-Planckian excursions of the inflation 
field for r < 2 X 10“® [12]. 

The first consequence of observing no tensors is there¬ 
fore that the model of inflation will be characterized by a 
potential of the form V (0) = Vb -f • • • with the constant 
first term dominating. 

The second consequence is that one can simplify the 
formulae for the spectral index and its running [2] 


Us-1 = 277*, as=2C*, (2) 

where r;* = MpV"jV and ^* = Mp{V'V'"jV'^). There¬ 
fore, observing no tensors implies that the running of the 
spectral index will be sizable only if the third derivative of 
the potential is nonvanishing. Also, since (rig — 1) = 2 ? 7 * 
is measured to be negative by Planck, we conclude that 
the second derivative of the potential must be negative. 

A prototype of inflationary models summarizing all 
these properties (domination of the vacuum energy, a 
nonvanishing H'" and V" < 0) is represented by the fol¬ 
lowing form of the potential (c > 0) 


H(0) = Vb (1 — c0^), —oo < p < 0 or p > 1. (3) 


This parametrization has the virtue to reproduce for 
p —>■ —oo the exponential models V = Vq (l — , 

among which the popular Starobinsky 7^^-model [13] 
and the Higgs inflation model [14]. Related to the ex¬ 
pression (3), as far as the predictions are concerned, 
are the logarithmic supersymmetry-inspired models V = 
Vb(l -I- c(ln0/p)) obtained for p —>• 0, and the brane in¬ 
flation models obtained for p = —4. Potentials of the 
form H(0) = Vo (l ± c0^) give a scale-independent spec¬ 
tral index. In all these cases the predictions are [2] (for 

P^2) 
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The running can be written as 
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where in the second passage we have made use of the 
known relation /nl = —5(ns — 1)/12 valid for single- 
field models [15] and the numerical estimate has been 
done using the Planck central value for (n* — 1) (as we 
will do in the following). The factor (p — 2)/(p — 1) is of 
order unity and the most favorable case is achieved for the 
logarithmic models, for which p ~ 0 and as — 1.2 x 10“®. 

Since the SPHEREx measurements of the power spec¬ 
trum and the bispectrum may provide Icr level detection 
of the running spectral index at the level of 6.5 x 10“^ 
[10], a non-detection of the running will rule out models 
with p < —12 at the Icr level. On the other hand a 2 ct 
detection of the running would require a running larger 
than 1.3 X 10“®, which is almost impossible to get in 
single-field slow-roll models of inflation if no tensors are 
observed, unless (rig — 1) is redder than the Planck central 
value by more than two Planck standard deviations. Of 
course, a large running can be obtained by abandoning 
the slow-roll regime [16], but in this case no correlation 
with the non-Gaussianity is expected. 

Multi-field inflation .— In the single-field inflation sce¬ 
nario the total curvature perturbation is generated dur¬ 
ing inflation and remains constant on super-Hubble 
scales. However, in many inflationary models there are 
usually a plethora of scalar fields which may play a sig¬ 
nificant role during inflation. Indeed, if these fields do 
not dominate the energy density during inflation, but are 
light enough to be quantum-mechanically excited during 
the de Sitter stage, they provide a source of isocurvature 
perturbation. 

If the curvature perturbation associated to the field 
driving inflation is suppressed, the total curvature per¬ 
turbation may be originated from such isocurvature per¬ 
turbation. From Eq. (1) one deduces that this is natu¬ 
rally achieved if the Hubble rate is small: no significant 
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amount of tensor modes are expected in the case in which 
the seeds of the LSS and CMB anisotropies are due to 
a light field other than the inflaton (this is in fact not a 
no-go theorem [17], but it is a very well-educated guess). 

We will concentrate on the most studied example of 
such a mechanism, the so-called curvaton mechanism 
[18-21], even though our findings hold as well for other 
mechanisms [22-25]. In the curvaton mechanism, during 
inflation the curvaton energy density is negligible and 
isocurvature perturbations with a nearly flat spectrum 
are produced in the curvaton field cr. 

After the end of inflation, the curvaton field oscillates 
during some radiation-dominated era, causing its energy 
density to grow and thereby converting the initial isocur¬ 
vature into curvature perturbation. The comoving cur¬ 
vature perturbation C can be computed through the dN- 
formalism [26] C{x,t) = N{x,t) — N{t), where N{x,t) 
is the amount of expansion along the worldline of a co¬ 
moving observer from a spatially flat slice at some initial 
time to a generic slice at time t. Since the number of 
e-folds from the end of inflation to the beginning of the 
oscillations is unperturbed because the radiation energy 
density prad = {ptot — Pa) dominates during this time, for 
the curvaton case one can redefine N as the number of 
e-folds from the beginning of the sinusoidal oscillations 
(that is when the quadratic part of the curvaton potential 
with mass m starts dominating) to the curvaton decay 
[27] 


N [pitd), p{tosc), o-*] = T In 
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We have highlighted the role of the function g = gicr^,). 
It represents the initial amplitude of curvaton oscillations 
which is some function of the field value at the Hubble 
exit. The generic function (/(cr*) parametrizes the non¬ 
linear evolution on large scales if the curvaton potential 
deviates from a purely quadratic potential with mass m 
away from its minimum [28, 29]. 

The resulting power spectrum of the scalar perturba¬ 
tions reads [27] 


Vs{k) 




( 8 ) 


where 


_ fg) 

ipaitd) +4prad(id) 

parametrizes the curvaton energy density at the time of 
curvaton decay with respect to the radiation component. 


We neglect for the time being the scale-dependence of 
77* and write d/dlnfc = Tfy^d/dt* = ((T*/i7*)9/9(j*. 
We then find 
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with r(| = (2p'/3(/)rd(l —rd)(3-|-rd). The expressions for 
the quadratic and cubic non-Gaussianities are [30, 31] 
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Notice that rd depends on cr* through Pa{td) in Eq. (7). 
This dependence is essential to obtain the correct /nl 
and pnl given above [27]. 

We can recast the running of the spectral index in 
rather appealing and compact form 
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The same expression can be derived more directly in 
terms of derivatives of the number of e-folds using the 
(5A7-formalism [26], but we have preferred to adopt a 
slightly more involved derivation to stress the importance 
of the non-linearities in the curvaton potential. 

Eq. (13) is the main result of this paper and indicates 
that the running of the spectral index can be sizable if 
5 nl > /nl- 1'^ general, this is true in those models for 
which rd — 1 and g ”9 j 9 '^ — 1 (or rd < 1 and g"gjg'^ — 
— 1) when /nl becomes small. This can happen if the 
curvaton potential has non-linearities in it [32, 33]. As a 
rule of thumb, being 3{ns — 1)^/4 ~ 9 x 10“^, we deduce 
that SPHEREx could detect a running of the spectral 
index at 2cr if (/nl/Znl = 0{±2). Since SPHEREx can 
detect at Icr-level a local /nl larger than 0.2, the situa¬ 
tion seems optimistic. 

Consider for instance the case in which the curvaton 
potential has a self-interacting piece Act^, where A is a 
dimensionless coupling, besides the quadratic term. In 
the limit in which the non-linear term dominates over 



4 




FIG. 1: Left: The running as as a function of /nl for a range of curvaton models. Right: Non-linearity parameter ^nl as a 
function of /nl for the same models as in the left plot. 


the quadratic piece during inflation g((T*) ~ cjV'^ [32, 33]. is a function of a light field cr. The corresponding power 
One obtains spectrum and non-Gaussianities are [31] 


A 2 c: 

ns-l = -A(rd(2-hrd)-l)^, /nl = ^ “ ’’d “ 2?’d) . 

(14) 

For instance, for Cd = 1/3 we get a red tilt, /nl = 5/9, 
(7nl — ~32/^l! Qfs ~ —13 X 10“^. This running can 
certainly be measured by EUCLID and SPHEREx. 

To describe further the correlation between a siz¬ 
able running and the level of non-Gaussianity, we 
present some illustrative examples obtained numerically 
in Fig. 1. For completeness, we have restored the scale- 
dependence of through an inflaton potential of the 
form U((/) ~ exp(/i(//Mp), which provides a constant 
contribution (—/i^/Mp) to (rig — 1). The various points 
have been obtained finding the function g{cr*)-, scanning 
over the various parameters and using the general expres¬ 
sion (11). We have systematically checked that all the 
observables, such as the amplitude of the power spec¬ 
trum, are in agreement with the current observations 
and that the contribution to the running from the V"- 
term in Eq. (13) and from the scale-dependence of iL* is 
smaller than the one coming from the ratio 3 nl//nl- 
2(7 SPHEREx forecasted errors for as and /nl are also 
shown. On the right panel we present the corresponding 
values of ^nl as a function of /nl- The dash line illus¬ 
trates that the running is mainly due to the correlation 
with the non-Gaussianity. 

Finally, let us consider the case of the modulated re¬ 
heating scenario. The curvature perturbation is gener¬ 
ated during reheating if the inflaton decay rate F = r(cr) 
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One can easily show that the analytical expression for 
the running is exactly the same as in Eq. (13). Now, at 
the time of reheating the held a will have a value cTreh = 
o’reh(o’*) and the plausible dependence of the decay rate is 
r(CTreh) cr/gjj {q > 0). If the potential for the light held 
is quadratic, then areh = c* and we find 5 nl//nl = 
the first term in the running turns out to be of the order 
of 6 X 10“^, too small to be detected. Nevertheless, if 
there is a non-linear relation between (Jreh and ct* due to 
the fact that the potential of the held a deviates from a 
purely quadratic potential, then one can easily increase 
the ratio 5 nl//nl and make the running detectable. 

Conclusions .—In this paper we have stressed that a siz¬ 
able running of the spectral index is typically connected 
to the non-Gaussianity in the case in which the curva¬ 
ture perturbation is due to the presence of more than one 
scalar field during inflation. A future detection of a signif¬ 
icant running of the spectral index and non-Gaussianity 
may tell us a lot about the true mechanism giving rise to 
the primordial perturbations. Dedicated searches aim¬ 
ing at measuring with great accuracy the running and 
non-Gaussianity are of vital importance. 
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